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The ionization energy theory is used to calculate the evolution of the resistivity and specific 
heat curves with respect to different doping elements in the recently discovered superconducting 
Pnictide materials. Electron-conduction mechanism in the Pnictides above the structural transition 
temperature is explained unambiguously, which is also consistent with other strongly correlated 
matter, such as cuprates, manganites, titanates and magnetic semiconductors. 
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The discovery of iron-based layered compound [l|,[l, 3| 
that has been found to be superconducting recently [4| 
have led to detailed electrical structural and mag- 
netic Q characterizations. The superconducting tran- 
sition temperatures (Tc) for this compound, namely, 
SmFcAsO0.85F0.15, LaFeAs0i_a:F2; and CcFcAsOi_a;Fa; 
arc in the vicinity of 42 K, 26 K and 30 to 43 K, re- 
spectively H 0, S B HO, [m, [II • One of the fundamen- 
tal requirements from the electronic characterizations of 
superconductors is to understand what parameters effec- 
tively control the flow of electrons, both above and below 
Tc. Similar to high- Tc cuprate superconductors [l^, the 
pairing mechanism in this new class of pnictide materials 
is also unknown presently [T^. Studies on this class of 
material may provide some additional clues on the over- 
all structure of the theory of superconductivity. When 
X is varied via substitutional doping, for example, sub- 
stituting F into O sites, two important but unexplained 
effects are observed, (i) systematic change in Tc and (ii) 
systematic change in the resistivity (or resistance) above 
Tc and above structural transition temperature (Tjv). 

In this work, we focus on point (ii) and evaluate the 
mechanism of electron-conduction above T/v for different 
substitutional doping, x in (Sm, Ce, La, Ba, Sr, Ca)- 
FeAsOi_:r-(F, CI, Br)^ and Bai_^-(Sr, Ca, K, Rb, Cs)^- 
Fe2As2 superconductors. Here, we provide unambiguous 
and consistent explanations on the electron-conduction 
mechanism based on the ionization energy theory 
which can be used to understand how and why both 
cationic and anionic substitutional doping change the re- 
sistivity and specific heat of these superconductors at 
temperatures above T/v. The advantage of this theory 
is that we do not need to know the microscopic details of 
the compound and therefore, the ionization energy the- 
ory (lET) is not suitable to evaluate the mechanism of 
superconductivity as stated in (i). 



Electrons are indistinguishable, but the energy of each 
electron is different in different solids and if we can ex- 
ploit this energy, then we can systematically study the 
changes of electronic properties in non free-electron met- 
als. This is also true for electrons in atoms due to dif- 
ferent magnitudes of electron-nucleus Coulomb force. In 
view of this, we propose here a methodology that inter- 
rogates the atoms to reveal information about their elec- 
tronic energies based on the many-body Hamiltonian, 



Hip - (Eo ± 0^. 



(1) 



H is the many-body Hamiltonian, while if is the many- 
body wavefunction. The eigenvalue, £'0 ± ^ is the total 
energy, in which Eq is the total energy at temperature (T) 
equals zero and ^ is the atomic energy-level difference. 
The -I- sign of ±^ is for the electron (0 +00) while the 
— sign is for the hole (—00 — *■ 0). Complete details about 
this theory and proofs can be found in Refs. [H, [3| and 
references therein. The average ionization energy value 
(after invoking the ionization energy approximation) can 
be obtained from 



Eli 
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The subscripts, i = 1, 2,...z, where z denotes the num- 
ber of valence electrons that can be excited or contributes 
to the electronic properties of a solid. In addition, the 
carrier density in) can be calculated from 
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where fe{Eo, Ej) is the ionization energy based Fermi- 
Dirac statistics, ks denotes the Boltzmann constant 
whereas, Ne{Eo) is the density of states. The specific 
heat formula is given by 
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where e is the electronic charge, A = (127reo/e^)as, 
h is Planck constant, c is the sound velocity, eo and as 
denote permittivity of free space and atomic hydrogen 
Bohr radius, respectively. Hence, all we need to know 
now is the relationship between Ei and x in order to 
predict the evolution of carrier density (conductivity) and 
the specific heat above Tm- 



III. ANALYSIS AND PREDICTIONS 

From the structural studies the oxidation (va- 
lence) states for the elements in (Sm,Ce,La)FeAsOi_i,Fa; 
are given by Sm^+, Ce^+, La^+, Fc^+, As^", and 
F~. Fluorine doped CeFeAsOi-^Fa; is equivalent to 
electron doping in which, the charge distribution be- 
tween the planes can be written as (FeAs)"^"* and 
((La,Ce)Oi_a;Fa;)^+'', where the 5 refers to the average 
number of electrons generated by Ce or La due to F 
doping. For example, 5 = \ for ce = 1, and 5 = if 
x = 0, where x can be related to any anions with 1 — 
valence state. These additional electrons then give rise 
to conduction in the FeAs planes d, Q. Our present 
understanding is that Cc or La will contribute to an ad- 
ditional electron in the (FeAs)^^^'' plane due to F dop- 
ing, which enhances the conductivity above T^f. How- 
ever, in order to discriminate the conductivity between 
Ce and La (or any other 3-1- valence state doping ele- 
ments), we need to employ the ionization energy the- 
ory introduced earlier. In addition, this theory can be 
used for any substitutional elements, be it anionic as in 
LaFeAsOi-a^Fa,, or cationic in both (Sm,Ce,La)FeAsO 
and Bai_:,K^Fe2As2 

Let us start examining the conduction mechanism in 
(Sm,Ce,La)FeAs0i_2^F^ compound (see Fig. [T|) due to 
increasing F (and decreasing O). Using Eqn. we 
have calculated the averaged ionization-encrgy values for 
all the elements discussed here arc listed in Table ID 
Prior to averaging, the ionization energies for all the ele- 
ments were obtained from Ref. [13] ■ We stress here that 
smaller Ei implies weak electron-phonon coupling that 
gives rise to easier electron-flow and also large carrier 
density p^ . However, this scenario is reversed if Ej is 
large. This means that systematic substitution of ele- 
ments with smaller Ei will systematically decrease the 
Ej of the compound and subsequently increases the car- 
rier density due to increased electronic excitation proba- 
bility. In principle, we can use this theory to fine-tune the 
conductivity of any strongly correlated matter. Firstly, 
we need to identify the source of the electrons, which is 
from the carrier doping layer, LaOi-^Fj,. Since the va- 
lence state for F is 1 — , and on average, this increases the 
electrons from LaOi-^F^ layer into the (FeAs)""'^"'' plane 
with respect to increasing F Q . The next question is 
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FIG. 1: Schematic crystal structure for 

(Sm,Ce,La)FeAsOi_a;Fi compound The source of 

the electrons (Sm, Ce, La) determines the ability to conduct 
electricity and these electrons availability comes from F 
doping. FeAs is the conducting layer, whereas (Sm, Ce, 
La)-(0,F) is the carrier doping layer. The 5e doping into the 
FeAs layer here is due to F substitution into O sites in the 
layer, (Sm,Ce,La)Oi-iFi,. 



TABLE I: Valence states and the ionization-energy values of 
selected elements (ordered with respect to increasing atomic 
number, Z). The average ionization energies were calculated 
from Eqn. ([2]). The unit kj/mol is adopted for numerical 
convenience. 



Elements 


Atomic number 
Z 


Valence 
states 


Ionization energy 
Ei (kJ/mol) 


F 


9 


1+ 


1681 


CI 


17 


1+ 


1251 


K 


19 


1+ 


419 


Ca 


20 


2+ 


868 


Fe 


26 


2+ 


1162 


Br 


35 


1+ 


1140 


Rb 


37 


1+ 


403 


Sr 


38 


2+ 


807 


Cs 


55 


l-f 


376 


Ba 


56 


2+ 


734 


La 


57 


3+ 


1152 


Ce 


58 


S+ 


1177 


Sm 


62 


3+ 


1292 



how does substitutional doping of CI or Br, instead of F 
will affect the resistivity and specific heat curves? Simi- 
larly, Sm or Ce substitution into La sites will also affect 
these curves systematically that needs to be understood. 
To answer these questions unequivocally, we need to in- 
voke the ionization energy theory, and its predictions are 
given in Fig. O 

We first discuss the electron conduction mech- 



3 



I (La, Ce. Sm)FeAs(0, F. CI, Br) 

(a 



Ida, Ce, SlTi)FeAs(0, F,Ci, Br|. 
(b) 




♦ (Ba Sr. Ca)Ff-,A.Si 

(c)' 




AiBa.K, Bb.C5lFi.:,.As.>- 

.(d) 



K Rb Cs 



FIG. 2: Predicted evolution of the resistivity (p) and inverse 
specific heat {Cv) curves (upward or downward) with respect 
to different doping elements (both cations and anions) based 
on the ionization energy theory. The predicted shifts depicted 
in (a) to (d) are valid for temperatures, T > Tn- The pre- 
dicted trends for resistivity and specific heat were obtained 
from Eqn. ([Sjl and Eqn. (|4]). Equation ((Sjl is used to average 
the ionization energies of all the elements. See text for details. 



anism and the evolution of specific heat in the 
(Sni,Ce,La)FeAs0i_2-F:i. and Bai_2:Ka;Fe2As2 com- 
pounds, for different doping elements. Since F, CI and 
Br are all anions, we are only required to know the first 
ionization energies of these atoms to evaluate the upward 
and/or downward shifts of the resistivity curve [Tsj . The 
reason that we used the first ionization energy is that 
the atoms with large Ej (F in this case) implies strong 
Coulomb attraction (due to smaller screening) between 
an electron and the ions. From Table [H we can write 
£:/(F) > £;/(Cl) > Ei{Bt), and after substituting this 
inequality into the ionization energy based Fermi-Dirac 
statistics [l5| and using Eqn. we can obtain the re- 
sults where the resistivity curve will shift downward (im- 
proved conductivity above Tn) with increasing content of 
CI or Br (see Fig. [2ja)). The reason is that when the el- 
ements with smaller ionization energies are substituted, 
then the electronic excitation probability will increase, 
which subsequently leads to increased carrier density or 
improved conductivity, and larger specific heat. There- 
fore, the resistivity curve (p(T > Tn)) shifts downward, 
while the specific heat curve shifts upward systematically 
with CI or Br doping as in 0i_2;Cl2; or Oi^x^t^x, as com- 
pared to Oi-xFx (see Figs.[2][a) and (b)). In other words, 
Br-doped samples will have the lowest resistivity, whereas 
F-doped sample will have the highest resistivity. The rea- 
son we used T > Tn is because when T < Tn, the elec- 
tronic excitation or the electronic polarization probabil- 
ity of the compound changes. Consequently, any change 
in these probabilities also corresponds to the change in 
the average valence states of all the elements in the com- 



pound. This is the reason why we see the structural 
changes (changes in the lattice parameters). 

On the other hand, if we were to substitute La with Ce 
or Sm, the resistivity and specific heat curves will shift 
upward and downward, respectively. In all the plots de- 
picted in Fig. [21 the y— axis denotes p{T > Tn), l/C-u 
where p cx 1/n oc 1/C„ based on Eqns. ^ and It is 
very important to realize here that there are different sets 
of Eqns. ^ and of which, one set of equations cap- 
tures the changes from F to CI doping in Fig.[2)^a) and an- 
other set of equations are required to capture the changes 
from CI to Br. The reason for this is that the electronic 
and phononic properties of a given sample in the presence 
of O and F has its own exponential trend that are differ- 
ent from the samples with O and CI, or O and Br. and 
so on for other elements shown in Figs. [2l^b) — (d). Add 
to that, in real sample preparations for different doping 
levels (a;), the valence states of all the elements may not 
remain the same for different x due to defects (interstitial, 
substitution at different sites and vacancies) [l6j . which 
further affect the exponential term given in Eqns. ([5]) 
and ([1]). As a consequence, all the solid lines connecting 
each calculated points in Figs. [5] and [3] have their own ex- 
ponential terms based on Eqns. ([3|) and (|4|). and therefore 
they are discontinuous. 

Using the same line of reasoning based on the ion- 
ization energy theory discussed above, we can also ex- 
plain and predict the effect of substitutional doping in 
Bai_a;(Sr,Ca,K,Rb,Cs)3;Fe2As2 compound. The upward 
shift in the resistivity and downward shift for the specific 
heat curves are clearly shown in Fig.[2jc) when Ba is sub- 
stitutionally doped with Sr or Ca. Moreover, Fig. [IJd) 
also indicates the effect of cation-doping at Ba sites, but 
these cations have 1-1- valence state in which, K, Rb or 
Cs will give rise to improved conductivity above Tn, and 
larger heat capacity, as compared to Ba. 
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FIG. 3: Theoretically calculated resistivities in (a) c— axis 
(shifts downward) and (b) afo— plane (shifts upward) for 
Lai_3;(Ca,Sr,Ba)a;FeAsO compound. Note that these predic- 
tions are observable if pj^ter > Pwri otherwise, pc oc pab- 



We stress here that the theoretical results presented in 
Fig. [2] are valid for both a6— plane and c— axis resistivi- 
ties, which implies that they are observable even for poly- 
crystalline samples. However, the analysis will become 
complicated for the Lai_a;(Ca,Sr,Ba)2;FeAsO compound 
due to different conductivities in FeAs and LaO layers. 
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In contrast, any substitution into Ba sites in BaFc2As2 
compound will affect both the a5— plane and c— axis re- 
sistivities equally (upward or downward: pab pc) be- 
cause there are no anions surrounding the Ba. Unlike 
BaFe2As2 compound, the LaFeAsO compound can be 
substitutionally doped in such a way to give contradict- 
ing resistivity curve shifts or trends between afe— plane 
and c— axis resistivities. These contradicting trends can 
be observed (i) if La is substitutionally doped with el- 
ements that have average valence states lower than 3-I-, 

(ii) elements with lower average ionization energies than 
La^+ are substituted and (iii) yof^y° > P^ytr- If condition 

(iii) is not satisfied, then pc oc pab, where both resis- 
tivities will have the same upward or downward shifts. 
For example, substituting La'^"'" with Ba^^ or Sr^+ or 
Ca^"*" in LaFeAsO compound will improve the electron 
conduction in the Lai_a;(Ba,Sr,Ca)a;0 layers, while the 
resistance in the FeAs layers will increase due to the lack 
of electron-doping from the Lai_a:(Ba,Sr,Ca)2;0 layers. 
Hence, by measuring the resistivities in afe— plane and 
c— axis of a single crystal, we will see a downward shift 
in the c— axis resistivity due to improved conductivity in 
Lai_3;(Ba,Sr,Ca)a;0 layers. Whereas, an upward shift is 
expected for a5— plane resistivities due to smaller elec- 
tron doping from Lai_a:(Ba,Sr,Ca)a;0 layers. The pre- 
dicted results based on the ionization energy theory are 
shown in Figs.[3](a) and (b). 

Thus far, there are two experimental results found to 
be in accordance with our predictions. Firstly, the resis- 
tivity curves of the polycrystalline Bai_:i-K3;Fe2As2 sam- 
ples indeed have a downward shift with doping x p^ . 
as correctly predicted in Fig. [l^d) (where K doped sam- 
ples have a lower resistivity at T > T^r). Apart from 
that, the polycrystalline Sr doped Lai_:ESr^FeAsO sam- 
ples shows an upward systematic-shift in the resistivity 
curve from x = 0.13 to 0.20 [2l|. This latter results indi- 
cate that Lai_i:Sri:FeAsO samples satisfy condition (iii), 
where p\^^j. > pf^J^f-- Based on the principle of least ac- 



tion, the resistivity measurements in polycrystalline sam- 
ples will be determined by the lowest resistance path [l^ , 
which corresponds to pf^^y in this case. Therefore, Sr 
doped Lai_2:Sr2:FeAsO single crystals are predicted to 
obey the theoretical results displayed in Figs, ^a,) and 
(b). 



IV. CONCLUSIONS 

In conclusion, we have provided a physical picture 
based on the ionization energy to predict and explain 
the electron conduction in (Sm, Ce, La, Ba, Sr, Ca)- 
FeAsOi_:i,-(F, CI, Br)^ and Bai_^-(K, Rb, Cs)^-Fe2As2 
superconductors above the structural transition temper- 
ature. Our predictions can be tested with the present ex- 
perimental techniques. Accurate predictions on the resis- 
tivity and specific heat has been made by identifying the 
type of ions (to obtain their respective energy- level differ- 
ence) and their valence states (to transform the atomic 
energy-level difference to the ionization energy) in the 
compound. By identifying the type of ions that exist in 
a given compound, we can evaluate the conductivity and 
the specific heat of the compounds accurately by means 
of the ionization energies. However, lET, in its present 
form, cannot be used to understand the origin of super- 
conductivity. 
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